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Abstract 


Electron-hole plasma heating and ultrafast modulation in a semiconductor 
laser under a terahertz electrical field are investigated using a set of hydrody- 
namic equations derived from the semiconductor Bloch equations. The self- 
consistent treatment of lasing and heating processes leads to the prediction of 
a strong saturation and degradation of modulation depth even at moderate 
terahertz field intensity. This saturation places a severe limit to bandwidth 
achievable with such scheme in ultrafast modulation. Strategies for increasing 
modulation depth are discussed. 
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I. INTRODUCTION 


As demand keeps increasing for faster data transfer rate in today’s broadband telecom- 
munication systems, conventional current injection modulation scheme on semiconductor 
diode laser and passive modulation schemes, such as those using electro-optic and acoustic- 
optic effects, are stretched to their limits . 1 As a result, both multiplexing technologies and 
alternative modulation schemes into the terahertz (THz) regime are under fervent pursuit. 
Direct optical modulation of the transmitter laser is distinctively advantageous over other 
passive schemes as the telecommunication industry is moving toward an all-optical network- 
ing future in order to eliminate the electronic bottleneck. One of such alternative direct 
modulation schemes is through a heating electrical field . 2 5 

In such a modulation scheme for quantum well lasers, a high frequency electrical field 
(HFF ) 2 ’ 3 or a terahertz field (THF ) 4,5 is applied to the laser structure such that the field 
polarization is in the plane of quantum well layers. The applied field accelerates electron- 
hole plasma and the acquired kinetic energy is then converted to thermal energy' through 
various scatterings to heat the plasma. For a time-dependent field, the plasma tempera- 
ture is expected to change, following the time-dependence of the applied field. Since it is 
well-known that optical gain is a sensitive function of the electron-hole plasma (EHP) tem- 
perature, the alternating plasma temperature will lead to an alternating optical gain, thus 
to a modulated laser output. In contrast to the conventional current injection modulation, 
where the bandwidth 6 , / 3dB , is limited by interband recombination carrier lifetime t c and 
density-induced differential gain dG/dN , modulation through a THz field should not expe- 
rience this speed bottleneck, since plasma heating is an intraband process which occur on 
picosecond time scale. The plasma density is almost constant on such a fast time scale, so 
that slow recombination is no longer an issue in this scheme. This is the basic idea behind 
the ultrafast laser modulation through a high frequency heating field. 

The basic phenomena of this modulation scheme have been discussed in several previous 
papers'’ 1 ’ ’ using simplified models. In those approaches, linearization of the optical gain at a 
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given temperature is made. The self-consistent coupling between busing and plasma heating 
processes is not fully considered. One obvious prediction of such a theory, though not 
explicitly emphasized in those papers, is that the modulation depth increases monotonously 
as a function of the THF/HFF amplitude. One is led to believe that the desired depth of 
modulation can always be achieved through an increase of the amplitude of the heating field. 
Since modulation depth is an important parameter for any practical implementation of this 
scheme, it is important to look at this issue more closely to find out if and to what degree 
the prediction of a simplified linear theory is valid. 

In this paper, we attempt to overcome these shortcomings of the linear theory and to 
address the two-way coupling between lasing and heating processes. While effects of plasma 
heating on lasing process has been addressed, the reverse process, namely the effect of las- 
ing on heating has not been adequately studied. Our starting point is the w’ell-established 
semiconductor Bloch equations (SBEs) that describe microscopically the coupling between 
plasma distribution and optical polarization induced by the laser field in the active semicon- 
ductors. We derive a set of equations for moments of up to second order from the SBEs after 
adiabatically eliminating the medium polarization. This set of moment equations consists of 
energy, momentum, and carrier density. Their self-consistent coupling and the microscopic 
gain are included automatically. 

This paper is organized as follows: In the next section, we will derive the set of mo- 
ment equations from the SBEs and review the basic assumptions we make in the process of 
derivation. In the third section, main results of our numerical simulation of these moment 
equations will be presented. Finally, the paper is concluded in the fourth section. 

II. HYDRODYNAMIC MODEL: THE MOMENT-EQUATION METHOD 

The physical system we are dealing with is a uniform EHP in a semiconductor quantum 
w'ell interacting with a single-mode laser field by stimulated emission and absorption, which 
is conveniently described by tin* two-band semiconductor Bloch equations . 7 9 The modified 
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SBEs 10 12 with an in-plane terahertz field in the direction of x, F T n : {t)=F T „,{t)x, arc used 
;us our starting point to derive the moment equations and listed below, 


0 <>F Tlli {t) 
dt h 

d _ e FTn z (t) 
dt h 

d tpTffz{t) 


n c (k) = -23 {%£-} + ^ 


Vf 


^ h 


(-«) = -23{n {P ;-}+*^il 


ih 


dt 


V £ ) + (&*>/ - £ e (k) - e h (k) - E s )] p * 
= hQg \n e {k) + n h (k) - l] + ih 


dt 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


where n C| y,(±£) is the nonequilibrium distribution function for electrons and holes with 
wavevector ±k, respectively. The sign S means the imaginary part of the complex quantity 
within the ensuing braces and the superscript symbol * means complex conjugate. The Q* 
is the renormalized Rabi frequency and defined through 


= d*(i)A + Y. v H'-kPk' (2.4) 

ic'^ic 

Denoted by A and p % are the slowly-varying envelop functions of the laser field, 
E(t)=Aexp(iuit)+c.c., and interband optical polarization, P^(t)=p^ exp(iuit) . For brevity, 
most notation of the time dependence has been dropped. As usual, e is the absolute value 
of electron charge, h is the Planck’s constant, Vg is the Fourier component of the screened 
Coulomb interaction between an electron and a hole, d cv (k ) is the interband optical dipole 
matrix element between an electron state, jcfc), and a hole state, j v-k), £ e ,h(k) is the carrier 
energy renormalized by the screened Coulomb interaction, huj t is the photon energy of the 
laser field, and E (J is the bandgap energy of the semiconductor gain medium. Major incoher- 
ent scattering processes, such as carrier-carrier and carrier-LO phonon collisions, are lumped 
in the terms labeled by subscript inc and treated in the second-order Born approximation. 13,0 
Last but not least, other contributions not included here will be added empirically to the 
moment equations, as appeared in Eqs. (2. 7-2. 9). 8 Before delving into further development 
of the above equations, it is helpful to discuss our treatment, of the scattering processes. In 
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this paper, temporal correlations are assumed to bo instantaneous for the statistical evolu- 
tion of the EHP. In the semiclassical picture, by neglecting memory effect, the incoherent 


scattering terms in the modified SBEs are expressed in terms of local carrier distributions, 
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where a=e,h, shorthands for electrons and holes, S is the active area of the QW. The 
Boltzmann scattering matrices \\ a (k,k') for electrons and holes and the scattering matri- 
ces W p {k, k') take into account of carrier-carrier and carrier-phonon scattering processes 
additively (see Ref. 14 for detail). 


It is computationally too expensive and very often unnecessary to solve the above k- 
resolved equations to describe a laser. Rather a set of moment equations is simpler and 
often used. This has been done for lasers without the heating field 15-18 , where the equation 
for the first-order moment, the momentum, was not necessary. Formally it is straightforward 
to write down equations for the various moments. For the first three moments, we have the 
following equations (see Ref. 19 for a general approach): 
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P„, E„) = §£ jl.A*.. } ”.„(k) , (2.10) 

ir is the quantum well width, and m r =m,,m/j/(m r 4- ni h ) is the reduced mass. The y- 
component of the momentum P n has no dynamical relevance in our case and thus been 
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left, out so P n stands for the r-componont. The upper/lower sign in Eqs. (2. 8, 2. 9) is for 
electrons and holes, respectively. Uniformity, or charge neutrality, of the system requires 
N,.=N h =N. Jinj is injection current with quantum efficiency rj. Details in carrier capture 
are not considered, but we assume that the capture process occurs instantaneously and pair- 
wise for electrons and holes. Energy gain per pair is intuitively taken as the bandgap offset 
■^Eg between barrier and well material, in addition to the average energy E a /N per pair 
injected into the EHP. 


It is worth noting that other than the introduction of those empirical terms accounting 
for carrier generation and recombination, no assumption has been made about the distri- 
bution functions up to this point. The collision terms, d{P a , E a }/dt\ Q0 , for electron-hole 
scattering, and d{P a , E a }/dt \ c _ LO , for carrier-LO phonon scattering, are treated in the Born 
approximation. 13 - 9 Their corresponding expressions can be easily obtained by multiplying 
hk x or h k 2 /(2m a ) to the scattering terms and then integrating over k. For example, 
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It is important to point out that the above set of equations is not closed yet. First, 
the distribution functions appear in the collision terms, Eqs. (2. 5, 2. 6). Assumption has to 
be made about the type and nature of these functions so that' those collision terms can be 
evaluated. Second, the ^-resolved polarization functions appear in the moment equations, 
Eqs. (2. 7-2. 9), which is an additional coupling that is not found in the typical moment 
equation treatment of electronic transport problems. This is the well-recognized truncation 


problem. 

Our way of truncating this hierarchy relies on the time-scale argument for the problem 
we deal with here. We note that the fastest time scale is the carrier-carrier scatterings 
which complete around 50~100 fs at the level of lasing densities. The next fastest time 


scah- is that of the carrier-phonon scatterings, which is in the order of a few hundred fs. For 
the terahertz field we consider here, the period is larger that the phonon scattering time. 
Since our main interest in this paper is on the terahertz heating and much slower lasing 
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processes, wo can safely assume that the carrier-carrier scattering processes have resulted in 
quasi-steady distributions for electrons and holes, or quasiequilibrium. This leads in general 
to the so-called drifted Fermi-Dirae distribution: 


«,(*) = /L°\k-k<g’) = {exp [p>M£- *<„">) - + l} _i (2,12) 


with hk = hk^px being the drift momentum of electrons/holes due to the applied 
THF/HHF in the x direction. And (3 = (kgT)~ l , Ub is the Boltzmann constant, T is 
the plasma temperature, c Q (£) is the unrenormalized carrier energy, and p a is the chemi- 
cal potential which is determined independently for electrons and holes from their density 
and temperature. 8 In doing so, the active semiconductor QVV medium is characterized by 
plasma density, temperature, and the drift momenta of electrons and holes, respectively. 
Here we have assumed that the plasma is thermalized by the subpicosecond incoherent 
scatterings such that no temperature difference exists between electrons and holes (see the 
Appendix). The assumption about the distribution functions cut off the hierarchy arising 
from the scattering terms. The same time scale argument also helps cut off the hierarchy 
due to polarization, since distribution thermalization and polarization relaxation are the 
results of the same carrier scatterings. This leads to a typical adiabatic elimination of the 
polarization. After defining p £ = e 0 x*:A and making a relaxation time approximation for 
the interband polarization scattering term in Eq. (2.3), 
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the fc-resolved optical susceptibility is written as 
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whore the term 1/(1 — Q\(k)) is the so-called Coulomb enhancement factor and its explicit 
expression can be found in Hamj and Koch* or in Chow and Koch 20 for the diagonal decay, or 
an improved version in Nincj cl. a l. {:} including the non-diagonal scatterings. The relaxation 
rate is assumed to be ^-independent. The many-body effects is treated in the? framework 
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of the Hat roo- Fork theory with dynamic screening modeled by the single plasmon pole 
approximation. 7 ' 8 It is worth noting the followings: (l) As discussed in the Introduction, 
subpicosecond incoherent scatterings allow us to neglect the coherent effects of the THz field 


on the polarization; (2) drift-related correction to the susceptibility has been omitted for the 
drift kinetic energy is estimated two orders of magnitude smaller than the chemical potential 
at a typical THF of a few KV/cm; (3) No spectral hole burning introduced saturation effects 
are considered in the optical susceptibility because the laser field intensity involved in our 
simulation is much smaller than the required intensity, estimated at 58 MW/cm 2 . 21 

If we substitute Eqs. (2.14) and (2.12) into the moment equations, we have a closed set of 
equations for the first three moments. This set, coupled with Maxwell equation, will permit 
a self-consistent description of lasing and heating due to the terahertz field. However, further 
simplification can be made to eliminate the momentum equations. First, we point out that 
in Eq. (2.8) the momentum change due to stimulated emission can be ignored in comparison 
with those due to terahertz field and to carrier-carrier and carrier-phonon scatterings. The 
momentum equation is then reduced to 

dP / p p \ 
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where 7 t Q is the lumped carrier momentum relaxation rate due to LO phonon scatterings, 

spontaneous emissions, and non-radiative recombinations, 7 e /, is the interband electron-hole 

scattering rate, which is given in the Appendix. Manipulation of this equation by Fourier 

transformation produces the following results, 
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Since the THz frequency u> is an order of magnitude smaller than j t +7 eft , 22 " 25 the frequency 
dependence is omitted in our investigation. 


The final procedure in arriving at the hydrodynamic model is to convert the energy 
equation, Eq. (2.9), to a temperature equation, which is aided by first decomposing energy 
E n into kinetic energy which is associated with drift momentum P a and thermal energy IF„. 
by using the following relation, 
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Then, the desired moment equations for the hydrodynamic model for plasma density N 
(zeroth-order moment) and total thermal energy (second-order moment), cou- 

pled to the laser-field amplitude-envelop-function equation, 26 are derived as follows, 
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where £ = \A\, K = y/eyui/c is the wavenumber of the single-mode laser field, c is the 
speed of light in vacuum, and V=wxS is the volume of the active region. In addition, 
t thz = 1/(t t + 7 e/i) is an effective relaxation time, similar to that in the Drude model 
for a single-component plasma. The weak energy dependence of both 7 1 and j eh has been 
neglected in our model, so a constant relaxation time has been used in this paper. Carrier-LO 
phonon interactions, — | , are treated microscopically and details are given elsewhere . 27 

Rather than solving the above set of equations for {£,N with the plasma temperature 
T as a dummy variable, it is more convenient to work with the temperature and numerically 
solve for {£,N,T}. A transformation is introduced below, which is equivalent to a similar 
manipulation , 16 based on the chain rule and the functional relationship, W=W(N, T), 
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where d\\ /dN and d\\ /OT are calculated using the assumed quasi-equilibrium Fermi-Dirac 
distribution functions for electrons and holes. 

III. NUMERICAL RESULTS AND DISCUSSIONS 


In this contribution, an S nin GaAs/Alo ..(Gao jAs single quantum well structure is stud- 
ied. Other than those material parameters, which can be easily found in literature and not 
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listed horn, we use the parameters in Table I. Throughout this study, we assume, except when 
specified otherwise, a fixed lattice temperature of 295 K and photon energy, huJi, of 1. 4 eV, 
which roughly corresponds to the gain peak photon energy. In addition, continuous-wave 
(C\V) THz field of the form, Frn z (t) = A' V uz sin (2Tns Tf [ z t — is assumed. The focus of 
this paper is on long-time behavior of the semiconductor laser under such CW terahertz 
field. The transient response and the effects of pulsing terahertz field will be considered in 
a forthcoming paper. 

As we pointed out earlier, the effects of THF/HFF on semiconductor laser modulation 
have been studied using a linear gain model. 3 ’ 5 According to the linear theory, a deeper 
modulation can always be achieved through further increase of THF/HHF amplitude. Since 
certain level of depth is required in application, it is important to know at what den- 
sity/temperature level, the linear gain model will become invalid. 

First of all, let us consider the nonlinear effect of application of an intense THz field on 
the differential gains induced by density change and temperature change. As understood 
from small-signal analysis in direct current-injection modulations, the important bandwidth 
parameter, f MB , is determined partly by the density-induced differential gain. When the 
optical gain modulation is realized via plasma heating, it is found that the temperature- 
induced differential gain replaces its density counterpart. However, deviation from linearized 
model is expected. The optical gain G of the Inverted medium, i.e ., the EHP, is given as 
below, 

c = (3-D 

In Fig. 1 (a) and (b), w r e plot the density- and temperature-induced differential gains 
{ dG/dN,dG/dT }, respectively. It is obvious from the figure that these differential gains 
change appreciably within the relevant density and temperature ranges (see Fig. 2 and 3). 
Such gain nonlinearities can be understood as follows. Because phase space filling effect of 
the fermionic nature of electrons and holes, saturation in the opt ical gain sets in as carrier 
density increases, which explains why the density-induced differential gain OG/ON is reduced 
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to zero. Plasma heating, however, reduces the phase filling factor, /,,(£) + /,,(£) - L, and thus 
eases the saturation and leads to a slight increase in dC/dN. On the contrary, as plasma 
heating reduces the phase filling factor, the temperature-induced differential gain OG/dT, 
proportional to the factor, decreases accordingly. As carrier density increases, saturation in 
the optical gain due to phase space filling suppresses the differential gain too. 

Now let us take a look at the asymptotic (long-time) behavior of the system under a CW 
THF. According to our extensive simulation in a quite wide range of modulation parameters, 
the asymptotic behavior of laser intensity, the EHP density and temperature, and optical 
gain are all well characterized as the sum of a zero-frequency (DC) component and a second 
harmonic (2i/ THz AC) component. Time evolution of these quantities start to deviate from 
this simple sinusoidal behavior when DC level of the laser output approaches zero — laser 
driven close to switch-off, or when the modulation frequency becomes too low and nearly 
resonant with the relaxation oscillation frequency. 5 This sinusoidal behavior with the second 
harmonic oscillation has its origin in the F% Hz {t) term in Eq. (2.20). The time evolution 
of laser intensity and carrier density can be found in Ref. 5. To look more closely at the 
modulation behavior as THF intensity changes, we plot the DC components (Fig. 2) and 
AC components (Fig. 3) of EHP temperature, EHP density, and laser intensity. There are 
quite a few features in Figs. 2 and 3 worth noting. First, all the DC components change 
with THF intensity linearly. As THF intensity increases, the DC plasma temperature rises 
due to heating effect. This temperature rise causes reduction in the optical gain. To tend 
to compensate for the loss of the optical gain due to temperature rise and to maintain laser 
operation, the EHP density has to increase to increase gain. This results in a net gain loss 
which cause a linear decrease of the laser output power as can be seen in Fig 2 (c). Finally, 
at the high level of the THF of our simulation, the gain loss due to the rising temperature 
far exceeds the density-induced increase, the net gain becomes comparable to threshold 
level and laser output drops to almost zero. The second feature of Fig. 2 is the frequency 
insensitivity of the DC components. As is obvious from the figure, the DC components for 
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all four frequencies, 0 . 1 , 0.33, 0.66 and iTIIz fall oti the same straight linos. This is because 
the inverse of the characteristic dynamical time scale in the EHP is at least an order of 
magnitude larger than these frequencies, thus the frequency dependence of the effective 
relaxation time T Tffz becomes negligible, as can be seen in Eq. (2.16). 

In contrast to the behavior of the DC components, the AC part shows a sensitive fre- 
quency dependence, as shown in Fig. 3. The AC components decrease as modulation fre- 
quency increases. This is a quite important point. Since a minimum modulation depth is 
required for application, this result places an upper limit to the modulation frequency. The 
most important feature in the AC component curves is the decrease of the AC components 
of density and laser intensity at high field amplitude (Fig. 3 (b) and (c)). While AC temper- 
ature shows a linear increase with THF intensity in the entire simulation range shown in the 
figure , the AC density and laser intensity increase only up to a particular THF amplitude 
when it reaches the maxima. Thereafter, the AC density and laser intensity start to drop. 
With further increase of the THF, the modulation depth approaches zero. This is contrary 
to the prediction of linear theory. 

To explain this behavior, we go back again to the quantity of optical gain. The total gain 
variation is written as the sum of two terms: AG = AG N + AG T with AG N = dG/dN x AW 
and £G t = dG/dT x AT. Typically in our simulation results, the AC component of density 
is in the order of AN = 10 8 cm -2 (see Fig. 3 (b)), while that of temperature is in the order 
of ten degrees (see Fig. 3 (a)). Using the typical number of differential gain in Fig. 1, we 
estimate that AG ^ < 0.1 cm~ l and AGt ~ 1cm -1 . This means that the gain modulation 
due to temperature modulation is at least one order of magnitude larger that that due 
to density modulation. This is actually the reason why this modulation scheme can get 
around of the interband speed bottleneck, since no appreciable density change is involved. 

As a consequence, we ignore the gain change due to density modulation in the following 
discussions. We note that the depth of laser output modulation depends on the ability 
to modulate laser gain, ft is therefore reasonable to assume that laser modulation depth 
is proportional to the gain modulation. For a linear variation of gain with temperature, 


12 


OG/OT is constant. With the monotonous increase! of AC temperature component as shown 
in Fig. 3 (a), a monotonous increase in AC laser intensity should be expected. The quite 
nonlinear behavior shown in Fig. 3 (c) indicates that the linear picture is no longer correct. 

To look at this issue more quantitatively, we plot in Fig. 4, the temperature-induced 
gain change as a function of AC temperature component using the linear gain model and 
a nonlinear model up to second order in AT. Since AC temperature is proportional to the 
THF intensity (Fig. 3 (a)), the x-axis can be also understood as the THF intensity. The 
parameters are exactly the same as in Fig. 3 with the reference temperature being 295 K. We 
see that, with the inclusion of gain nonlinearity, the gain modulation first increase up to a 
maximum and then starts to decrease as the nonlinearity starts to kick in. This explains why 
the AC laser intensity decreases as seen in Fig. 3 (c). This result has important consequence 
for application of this modulation scheme. Instead of simply increasing the THF intensity 
to increase modulation depth, a more careful optimization is required to achieve maximum 
modulation depth. This is especially important, since the relative modulation depth in this 
scheme is quite small very often. An estimate from Figs. 2 and 3 shows that the relative 
modulation depth is around 10 -5 . 

Indeed, maximizing the modulation depth is an important task for the application of 
this modulation scheme. From our numerical results presented above, we can point out 
several directions to pursue for this maximization. One obvious direction is to decrease the 
modulation speed as the modulation depth increases inversely with frequency as can be seen 
in Fig. 3 (c). In fact, we have not been able to achieve a reasonable modulation depth 
for frequency beyond 1 THz with CW THF of a few KV/cm. With careful optimization, 
this scheme mav achieve decent depth of modulation in the range between a few tens of 
gigahertz (GHz) to around 100 GHz. This is still a very appealing bandwidth, given the 
fact that the direct current modulation can achieve around 10 GHz. Another direction for 
maximization is to optimize the laser operating point, i.<\, the operating threshold current 
density etc. Plotted in Fig. 5 are long-time simulated dynamics of the output intensity of 
a GaAs Q\V laser operating at 77 K and under higher injection levels, / and 2 KA/nn 2 . 
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It, is dearly indicated that modulation depth increases with the injection current. As a 
result, a decent ~ 5 % relative modulation depth can be achieved at high enough modulating 
THF intensity. It should be pointed out that the increase in modulation depth mainly 
depends on the injection level instead of the lattice temperature. From our analysis above, 
we understand that the temperature-induced gain nonlinearity is the major obstacle for 
achieving deep modulation. From Fig. 1 (b), we see that dG/dT becomes a less sensitive 
function of temperature (meaning a smaller gain nonlinearity) as density increases. Increase 
of lasing threshold will have the laser operate around larger carrier density and therefore 
higher bandwidth. This can be seen in Fig. 4 by comparing the two curves for N=l.O x 
10 12 cm~ 2 and iV=1.8 x 10 12 cra -2 . Much less severe saturation of gain modulation is seen 
for the larger density case. Therefore increase in modulation depth is expected for the 
larger density. There is another possibility to maximize the modulation depth. We know 
that the gain characteristics depend on material properties and quantum well structures. 
By engineering quantum well structures with smaller gain nonlinearity, we may be able to 
increase modulation depth. Finally, modulation depth can be increased if pulsed THz field 
is used. This and other related issues will be discussed in detail in a forthcoming paper. 

IV. CONCLUSIONS 

To conclude, we derived a set of moment equations from the microscopic semiconductor 
Bloch equations describing a semiconductor quantum well laser in an in-plane THz electric 
field. This self-consistent treatment of the lasing and heating processes in the presence of 
a heating field allows a better treatment of this system. As a result, we predict the strong 
saturation and degradation of modulation depth at relatively small THz intensity level. 
This places a severe limit to the bandwidth achievable using this scheme. The saturation 
phenomenon is fully understood using the nonlinear gain model, which we obtained auto- 
matically by starting from tin* microscopic theory. Based on our understanding, we are able 
to point out several strategies for achieving deeper modulation. While this heating scheme 


proves not to he viable for frequency beyond 1 THz, it provides an attractive alternative at 
100 GHz level or lower. 

The key issue with this scheme that has to be resolved is how to apply the THF or 
HFF to a semiconductor laser. There are two ways to realize this scheme. Adding lateral 
electrode to bias the laser structure in the lateral direction provides an easy way as was 
suggested in Ref. 3. For demonstration purpose, terahertz wave can be coupled from the 
side of the laser structure into the active region. The mathematical treatment presented in 
this paper should be applicable for both cases. 

Finally, we point out that semiconductor laser driven by a lateral heating field may 
have other application, such as switching, HFF detection, and short pulse generation etc. 
The detailed discussion of these issues will be beyond the scope of this paper and will be 
presented elsewhere. 
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APPENDIX: MOMENTUM RELAXATION DUE TO ELECTRON-HOLE 

SCATTERINGS 

Carrier-carrier collisions occur on the femtosecond time scale at typical lasing plasma 
density. 2 ' 1,25 Intraband and interband scattering rates are comparable. To investigate the 
consequences of interband scatterings on momentum relaxation we assume quasiequilibrium 
for the electron-hole plasma. 

Given that each subband is thermalized and described by a drifted Fermi-Dirac distri- 
bution function, = // t ()) ( k — kff'j, where a—eji, the electron-hole scattering rate can be 
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written as follows 
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- n 2 n U l ~ n l)( l ~ n 'LJ 


k k* +(p * ‘ ' v ‘ ' v k+q* 

where a and /3 differ, and 1} is the screened Coulomb potential. After factorization, the 
difference term in the squared bracket contains a factor 


A 


3 


(1 -exp[- 


k' +q~k^ D ^ k' -k.Q 


k-kl 


... . - hq-(!&- 

‘ + ’-‘l]) = (1 - cxp(-C^L 


e? 


m h 


]) 


k B T J/ v ~ ~” t ' 1 k B T 

where we have used energy conservation and parabolicity of carrier energy dispersion. Ex- 
pansion of the exponent term leads to the conclusion that 


dn\ 

k 


dt 
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hk e D 

hk h D _ 
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II 

1 

<** , 

cx 

eh 

m e 

m h 

U 

“mj 


X . 


(A3) 


After integrating over k for the scattering terms above according to Eq. (2.11), and 
realizing that the following relationship holds under the assumption of quasiequilibrium, 


P a = Nhk a D , 


(A4) 


it is straightforward to show that the momentum relaxation due to electron-hole scatterings 
can be written in a rate equation form as in Eq. (2.15) 


dP e 


dt 


dPh 


eh 


dt 


~7ehm r 


eh 


(Pe_ _ Ph_\ 

\m e m' h ) 


(A5) 


and the relaxation rate is given by 

2y 


7 eh — ^2 ^rrVffi[ £ e{k)+£h{k' + q)—£ e (k + q)—£h(k')] 

ti 0 ) (k)fl 0) (* + ?) [i - /i 0) (£ + $)] [i - /f(^)] 


(A6) 


r hk 


h k^ 

TTlfr 


in the leading order of m 

Finally, we should comment on the use of the equality T e = 7), in deriving Eq. (A3). 
Exactly speaking, T e = T \ should be reached on the same time scale as in Eq. (A5). This 
means that the above derivation is not self-consistent on time scale of . Since we finally 
eliminate P,. and P/, adiabat.ically using Eq. (2.1G), our treatuK'iit on long time scale is, 
however, self-consistent. 
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TABLES 


TABLE I. Partial Simulation Parameters 


r 

confinement factor 

0.025 

T ph 

photon lifetime 

40 ps 

Jinj 

injection current 

0.4 KA/cm 2 

V 

current inj. efficiency 

0.5 


spontaneous emission 

X 

Cl 

1 

0 

Tnr 

non-radiative recomb. 

I.OxICP 


coefficient 

cm? /s 


coefficient 

s~ l 

TTHz 

relaxation time 

OA ps 

hip 

dephasing rate 

5 meV J 
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APPENDIX: FIGURE CAPTIONS 


Fig. I: Differential gains versus temperature at various densities. Density-induced differ- 
ential gain and temperature-induced differential gain are plotted in (a) and (b), respectively. 

Fig. 2: Dependence of DC components of long-time solutions for plasma temperature 
(a), density (b), and laser intensity (c) on THz field intensity at 0.1, 0.33, 0.66, and 1.0 THz. 

Fig. 3. Dependence of modulation amplitudes {2vi component) of plasma temperature 
(a), density (b), and laser intensity (c) on THz field intensity at 0.1, 0.33, 0.66, and 1.0 THz. 

Fig. 4. Temperature-induced nonlinearity in optical gain. Plotted is the change in gain 
due to plasma heating, fcG — dG/dT\ Tr x ZXT + d 2 G/dT 2 \ Tr x (zXT) 2 / 2, where reference 
temperature T r is 295 K, at two plasma densities. Dashed-line data only keep the linear 
term. 

Fig. 5: Long-time behavior of laser output operating at 77 K and under THz field 
modulation. Injection current for the dotted curves is 2 KA/cm 2 and 1 KA/cm 2 for the 
solid curves. 
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